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The worldsheet formulation of lattice gauge theories has two appealing features: the gauge non-redundancy 
and the geometrical transparency. Both properties are profitable in order to perform numerical computations. In 
the case of dynamical fermions this description offers additional advantages. For instance, it does not suffer from 
the species doubling problem and it involves fewer degrees of freedom. 



1. The Worldsheet Formulation Approach 

Simulating QCD with dynamical fermions or 
full QCD is still too expensive in computer time. 
Year after year the conclusion at lattice meetings 
is the same: either more machine power and/or 
a real improvement in algorithms are needed to 
produce reliable estimates of hadron spectroscopy 
We propose to explore a new numerical ap- 
proach based on the socalled F '-representation 
Staggered fermions can be naturally introduced 
by considering open paths representing string- like 
"electromeson" configurations in addition to the 
closed ones or loops representing pure electric flux 
excitations. This gives rise to the space of paths 
P of the above mentioned P-representation 
All the physical operators have a geometrically 
transparent realization in this space. 

The worldsheet formulation of lattice gauge 
theories (in terms of variables attached to the 
worldsheets of P-paths) was introduced for sev- 
eral models ||-f3|. U(l) theory || and scalar 
electrodynamics (SQED C ) Q were simulated us- 
ing this representation and the results were quite 
encouraging. More recently, a worldsheet descrip- 
tion of gauge theories with dynamical fermions 
was proposed fi. 

The worldsheet partition function Zp of QED 
is expressed as a sum over the worldsheets of 
strings or paths of the P-representation. That 
is, surfaces such that: (I) their borders T c are 
self-avoiding polymer-like loops and (II) when in- 
tersected with a time t — constant plane they 
produce paths begining at even sites and end- 



ing at odd ones. This description differs from 
that one gets by integrating the fields in the 
Kogut-Susskind action in two features. In first 
place, the isolated links traversed in both oppo- 
site directions or "null" links were ruled out and 
with them the myriad of different configurations 
corresponding to a given configuration of world- 
sheets. In second place, by virtue of the con- 
straint (II), the surfaces when intersected with 
a time t = constant plane produce only paths 
begining at even sites and ending at odd ones 
instead of the paths of ends with whatever par- 
ity. This two differences translate respectively in: 
computer time economy and a cure of the species 
doubling problem of the Kogut-Susskind action. 
The expression of Zp is as follows M : 



Zp = ^>(.nexp{-^ £ n% 



(1) 



pes? 



♦Supported in part by CONICYT, Project No. 318. 



where the sign ajr is given in terms of purely ge- 
ometric quantities of the fermionic loops T c || . 
In ref. || it was proved that Zp leads to the QED 
Hamiltonian using the transfer matrix procedure. 

Our aim was to test the worldsheet formula- 
tion of ref. H as a method to simulate dynam- 
ical fermions. Hence, we chose the simplest lat- 
tice gauge theory with fermions, the Schwinger 
model or (1+1) QED, and we performed a Monte 
Carlo simulation. This massless model can be 
exactly solved in the continuum. However, it 
is rich enough to share many features with 4- 
dimensional QCD as the confinnement and the 
breaking of the chiral symmetry with an axial 
anomaly. For this reason it has been extensively 
used as a laboratory to analyze the previous phe- 
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nomena and also to test different techniques to 
simulate theories with fermions. 

2. The Schwinger Model in the Worldsheet 
formulation 

In D = 2, it turns out that <j(J- c ) can be ex- 
pressed in terms of the number of connected parts 
Njrc, the lenght Ljr c and the area Ayrc of T c 
as a{T c ) = (-l) N ^-^+ A ^ and all the non- 
vanishing contributions have a? = +1. The rea- 
son is that Njrc — — + Ajrc = the number of 
enclosed sites by the fermionic loops T c which is 
always even by virtue of the above constraints (I) 
and (II) ( see ref.Q for more details), so that we 
omit this factor. The fermionic paths T c can be 
expressed in terms of integer 1-forms -attached to 
the links- / with three possible values: and ± 1 
with the constraint that they are non self-crossing 
and closed (df = 0) as 
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The lattice chiral condensate per-lattice-site 
is defined as < \X >= 2N;^x(~ l ) Xl < 
[x^(x),x( x)] >, where N s is the number of lat- 
tice sites. The corresponding operator is realized 
in the P-representation and thus we get for the 
chiral condensate j|: 



< XX >= 2 ~ 



2Np 
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where Np is the number of connected paths at 
a given time t. Thus, equation (||) allows to cal- 
culate directly the chiral condensate simply by 
counting the number of "electromesons" we have 
when we intersect their world sheets with each 
time slice t. 

3. Monte Carlo Simulation and Results 

To generate the worldsheets of the string-like 
configurations we use a Metropolis- type updating 
algorithm with the Boltzman weight proportional 
to exp{ — J2 P n P }- We simulate the model on 
Lx L square lattices with periodic boundary con- 
ditions (pbc). 



All the plaquettes p belonging to an open sur- 
face, by virtue of the self-avoiding constraint on 
their borders, must have all their plaquettes ful- 
filling the condition that the difference between 
integers n p of contiguous plaquettes has 3 possi- 
bilities: An p = ±1 or 0. There are also closed sur- 
faces corresponding to the world sheets of loop- 
like pure electric flux lines. In two dimensions the 
only possibility for a pure electric (closed) surface 
is to cover the entire lattice and then being closed 
by the pbc. Thus, all the plaquettes p belong- 
ing to it have assigned the same integer value of 

n p — ±1, ±2, Our algorithm works as follows: 

whenever there are different possible values for 
the integer n p attached to plaquette p - such that 
the produced surface fulfills the constraints (I) 
and (II) (self avoiding frontiers and odd length 
string respectively) - the present value is changed 
to one of these values and the Metropolis test is 
performed. Sometimes there is only one possible 
value of n p consistent with the constraints (I) and 
(II), in such cases the n p is updated to it and 
the Metropolis test is skipped. Typically, 100000 
sweeps per point were performed. 

To compute < \X > using (||) we have to count 
Np at each one of the L time slices. Therefore for 
each lattice sweep we collect L values of < \X >■ 

First, we checked that we get the right strong 
coupling behavior both for the ground-state en- 
ergy density ojq and the chiral condensate per- 
lattice-site < XX >■ The agreement with the se- 
ries expansion is pretty well up to f3 — 0.12 as 
one can see from FIG.l. 

On the other hand, the continuum theory is 
reached at zero lattice coupling, in the same way 
as four dimensional asymptotically free theories 
like QCD. It is known exactly that the chiral con- 
densate in the continuum is given by 



< 4>ip > 



= 0.15995. 
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As long as the lattice size L is increased, the 
convergence of q =< XX >/e to its known con- 
tinuum value (Q) improves. A property of the 
action of (||) is that q becomes stabilized at its 
continuum value well inside the weak coupling. 
Before reaching this regime, this observable oscil- 
lates strongly with (3. In FIG. 2 we plot q vs. (3 
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Figure 1. (a) The strong coupling behavior of io Q 
on a 4x4 lattice (up triangles), (b) The corre- 
sponding strong coupling series (down triangles). 



Figure 2. q —< XX >/ e f° r L —16 (squares), L = 
20 (diamonds), L = 24 (triangles) and L = 28 
(circles). The horizontal line correspond to the 
continuum limit value ~ 0.16. 



for lattices of 16x16, 20x20, 24x24 and 28x28. 
The value at which the simulation stabilizes at 
the weak coupling regime decreases from around 

0. 29 for L = 16 to 0.17 for L = 28. 

4. Conclusions and Final Remarks 

Our proposal was to show that the recently in- 
troduced worldshcct formulation was a valuable 
and alternative tool in order to do numerical com- 
putations with dynamical fermions. The method 
present the following advantages: 1) The easi- 
ness of computation. 2) Economy I: no gauge 
redundancy. 3 Economy II: it involves fewer de- 
grees of freedom that the Kogut-Susskind action. 
4) As a consequence of the constraint (II), the 
worldsheet action does not suffer from the species 
doubling problem of the Kogut-Susskind action. 

The values for the chiral condensate are consis- 
tent with the exact results for the continuum the- 
ory (recovered at the weak coupling fixed point). 
In addition, our strong coupling results are in 
agreement with the strong coupling expansion. 
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